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ABSTRACT 



Assuming that the electroweak and QCD phase transitions are first order, upon 
supercoohng, bubbles of the new phase appear. These bubbles grow to macroscopic 
sizes compared to the natural scales associated with the Compton wavelengths of 
particle excitations. They propagate by burning the old phase into the new phase 
at the surface of the bubble. We study the hydrodynamic stability of the burning 
and find that for the velocities of interest for cosmology in the electroweak phase 
transition, the shape of the bubble wall is stable under hydrodynamic perturba- 
tions. Bubbles formed in the cosmological QCD phase transition are found to be 
a borderline case between stability and instability. 
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1 . Introduction 

There is renewed interest in the dynamics of the electroweak phase transi- 
tion [1-4]. This is due in part to the possibihty that the baryon asymmetry of 
the universe might be generated at such a transition [5-9]. If the phase transition 
is first order then all of the conditions necessary for generating a baryon excess are 
present in the electroweak theory and its generalizations. 

The Sakharov conditions for generating a baryon number excess are [10]: 

1. Baryon Number Violation 

2. CP Violation 

3. Lack of Thermodynamic Equilibrium 

The first condition is now well established in electroweak theory [11-16]. The 
rate of baryon number violation at temperatures larger than T > 100 GeV is known 
to be much greater than the expansion rate of the universe, and not much less than 
a typical particle scattering time. 

The second condition is satisfied in the standard model, although explicit cal- 
culations show that the effects associated with CP violation arising from the quark 
mass matrix is much too small to generate an acceptably large baryon asymmetry. 
This is straightforward to patch up in generalizations of the standard model by 
allowing CP violation in the scalar interactions in an extended Higgs sector. 

The third condition is satisfied if the electroweak phase transition is first order. 
The transition to the ordered low temperature phase takes place by bubble nucle- 
ation. Since all of the old phase must pass through the walls of a bubble where 
the matter is strongly out of thermal equilibrium, all of the matter is strongly out 
of equilibrium at some time during the phase transition. If the baryon asymmetry 
is made near or at the bubble walls, then all of the Sakharov conditions can be 
satisfied. 

The electroweak phase transition is believed to be of first order. This has been 
shown to be true to the first few orders in perturbation theory [2-4]. In higher 
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orders, if apy is small enough, then these calculations are reliable. For realistic 
values of ~ 1/30, these results are somewhat controversial, and perhaps may 
only be definitely resolved using non-perturbative lattice methods. In this paper 
we will assume that the electroweak phase transition is of first order, and will 
accept at face value the results of perturbative computations. 

If the QCD hadronization phase transition is of first order, there might be 
consequences in both cosmology and heavy ion collisions. Whether or not the 
transition is first order is the subject of much study using lattice non-perturbative 
methods [17]. The order of the transition depends sensitively on the number of 
quark flavors and the quark masses and the physical case seems to be at the 
borderline between first and second order transitions [18]. In the following, we 
shall assume that the QCD phase transition is first order. 

To proceed with a semi-quantitative analysis of the effects of a first order 
phase transition either in cosmology or during a heavy ion collision, one needs 
to understand the formation and growth of phase transition bubbles. For the 
electroweak and QCD cases, much is now known [19-25]. In particular, the limiting 
velocity of such walls has been computed for electroweak theory, and the shape of 
the wall near the burning surface is known. In the QCD case, the macroscopic 
features of the bubbles have been studied. 

In this paper, we discuss the hydrodynamic stability of the bubble propaga- 
tion. This has been discussed previously for the QCD phase transition [26]. We 
find however that the previous analysis was incomplete. If this analysis is applied 
to either QCD or electroweak theory, one would conclude that whenever the wall 
velocity is less than the sound velocity of the matter, then the bubble wall prop- 
agation is unstable. Our results are in contradiction with these conclusions. The 
source of the discrepancy is a more thorough treatment of the boundary conditions 
for the motion of the bubble wall. 

The outfine of this paper is the following. In the second section of this paper, 
we derive the macroscopic features of phase transition bubbles for the electroweak 
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phase in the 1 + 1 dimensional approximation, extending known results for the case 
of QCD [20-21]. In the third section, we derive the equations which describe a 
possible convective instability which could affect bubble wall propagation. Our 
analysis is fully relativistic, but we also discuss the non-relativistic limit. We point 
out the origin of the discrepancy between our results and those which have been 
previously derived. In the fourth section, we solve for the possible instability for the 
electroweak and QCD transitions. For the electroweak theory, we find stable growth 
for all velocities larger than s/awj^ ~ 0.07. Detailed computations [23,24] of 
velocities favor v > 0.1, and so this instabihty is of no physical interest. For the 
quark-hadron phase transition, our analysis is inconclusive. A conclusive statement 
can be made only after the relevant parameters of the phase transition are known 
more accurately. In the last section, we summarize our results. 

2. Macroscopic Bubble Features 

Suppose we have a system which has a first order phase transition at a tem- 
perature Tc- If the system is expanding, then the system supercools below Tc 
before nucleating droplets of the low temperature phase. After these droplets have 
formed, they begin to grow and expand into the system. If the supercooling is suf- 
ficiently strong, as is the case in the electroweak phase transition, then the droplets 
expand until they collide with one another and complete the phase transition. If 
the supercooling is sufficiently weak, as is the case in the QCD phase transition 
in cosmology, then the droplets expand until there is a mixed phase of droplets 
of the low temperature phase embedded within the high temperature phase. The 
subsequent evolution of the mixed phase involves coalescence of droplets, and as 
the system expands, an increase in the volume of the low temperature phase at 
a rate determined only by the rate of expansion of the system. Eventually, all of 
the system falls into the low temperature phase when expansion has reduced the 
energy density of the system to that of the low temperature phase. 

In either case, for some time after nucleation, the bubbles freely expand into 
the system eating up the high temperature phase. If the growth of the bubbles 
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is hydrodynamically stable, then we expect that the bubble shape is smooth over 
macroscopic distances and should be spherical. The microscopic physics is only 
important in a region of order the diffusion length, a typical particle mean free path, 
near the surface of the bubble wall. This is seen explicitly in the hydrodynamic 
equations, which admit similarity solutions, that is, where the typical size scale is 
proportional to the time [20] . 

The analysis of a realistic situation is however a little complicated. There are 
two cases one must distinguish, referred to as deflagrations and detonations. 

A detonation occurs when the velocity of the burning front is larger than the 
sound velocity of the matter in the high temperature phase. This turns out to be 
equivalent to the case when in the rest frame of the burning front, the velocity of 
matter behind the burning front is less than that in front of it. 

This burning front in its rest frame is shown in Fig. 1. This is a frame which 
we will use for analyzing possible convective instabilities. For understanding the 
physics, it is sometimes more useful to consider the reference frame where the 
bubble wall moves, but the matter far in front of it is at rest. 

For the single front shown in Fig. 2 we have four unknowns, all other thermo- 
dynamic quantities follow from the equation of state, p = p{T), s = p'{T), e(T) = 
Ts{T) — p{T). We shall neglect all chemical potentials. Between the four un- 
knowns, two equations follow from conservation of energy and momentum fluxes. 
By requiring that d^Q^'^ — with* Q^'^ — + p)u^u'^ — pg^^ , it follows that 
A©^^ — 0, where z is the direction of motion of the bubble wall. From these one 
can, for instance, solve for the two velocities in terms of the two temperatures as 



The system is then completely determined by firstly giving the temperature Tq 
of the supercooled q-phase into which the bubble is expanding. This is obtained 

★ Here u** is the fluid four-velocity satisfying v? = 1. Our metric has signature (1, —1, —1,-1). 




(Ph -Pq)(eq +Ph) 

(eh-eq)(eh+Pq)' 



,2 ^ (Pq -Ph)(gh +Pq) 

" (eq-eh)(eq+Ph)' 



(1) 
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from a nucleation calculation. Secondly, a microscopic calculation [23-25] gives the 
velocity fb(T'q) of the combustion front. This has so far been calculated for the 
case of — -C . 

Boundary conditions often demand initial or final matter (or both) to be at 
rest and the single front discussed above is not sufficient. This leads one to a study 
of bubbles [20-21]. To understand a detonation bubble it is useful to introduce 
space-time rapidity, 



A 1-time -|- 1-space-dimensional detonation bubble is shown in Fig. 2, together 
with a plot of the flow rapidity and energy density as a function of space-time 
rapidity. Notice that constant space-time rapidity corresponds to x/t ~ constant, 
that is, to similarity growth. It is well known that the hydrodynamic equations 
admit similarity solutions. Since the front is traveling faster than sound velocity, 
the matter in ahead of the front is unaffected by the front until the front actually 
comes in contact with the matter. Behind the burning front is fast moving com- 
pressed matter, which over some region of rapidity has a constant flow velocity and 
energy density. The flow rapidity must be determined by solving the hydrodynamic 
equations. Behind this region of constant flow velocity is a similarity rarefaction 
wave at which the flow rapidity goes as Q{y) = y — ys, where ys — tanhcg is the 
rapidity associated with the sound velocity Cg. The flow rapidity can vanish when 
the space-time rapidity is the rapidity of sound. Therefore this point moves with 
the sound velocity and tries to catch up with the front. It can never do so, because 
the front is moving supersonically. The point at which the similarity rarefaction 
wave begins is when y — yi — ynow + Vs- In the very center of the bubble, the fluid 
velocity must be zero. 




(2) 



and flow rapidity 




(3) 
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A deflagration bubble is shown in Fig. 3. A deflagration (slow burning) bub- 
ble propagates into the high temperature phase at a velocity less than the sound 
velocity of the high temperature phase. This turns out to be equivalent, in the 
rest frame of the burning front, to having the velocity of the fluid behind the front 
larger than that ahead of it. In this case, the matter in front of the bubble is 
compressed and accelerated to a finite velocity. In front of the burning front, the 
moving matter produces a shock which moves faster than the sound velocity. (This 
is analogous to the shock front produced by moving a piston down a pipe) . Inside 
the burning front, the matter is at rest and at a constant density. Outside the 
shock front the matter is at rest and at constant density. The distance between 
the shock front and the burning front increases as a function of time. 

To compute the properties of these detonation and deflagration bubbles, one 
must solve the hydrodynamic equations. In the regions between but excluding 
the fronts, there are acceptable solutions of the hydrodynamic equations. For 

the deflagration bubbles, the region excluding the front is just a fluid moving at 
constant velocity, which will solve the hydrodynamic equations. For the detonation 
bubbles, the rarefaction region has a solution which we will later explicitly compute 
for the electroweak phase transition. 

Before explicitly solving the bubble equations, it is useful to review the proper- 
ties of the solutions. For deflagrations, we know the energy density and velocity of 
the fluid to the right of the shock front. We know the velocity of the fluid inside the 
burning front. We must determine the velocity of the burning front and the shock 
front, the energy density inside the burning front, and the velocity and energy 
density of the fluid in the region between the shock front and the burning front. 
These are 5 unknowns. The equations for the continuity of energy and momentum 
flux across each front gives four equations. To solve the problem we need one more 
equation which is microscopic. In electroweak theory, this is the equation for the 
Higgs fleld. 

For detonations, we know the velocity and energy density of the fluid to the 
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right of the burning front. The energy density and velocity of the fluid to the left 
of the burning front and the velocity of the burning front are the 3 unknowns. 
The continuity of energy and momentum flux are 2 equations, and again we need 
a microscopic equation. 

The bubble solutions have been explicitly discussed in Ref. 20 for the 1 + 1- 
dimensional case and in Ref. 21 for 1 + d = 1, 2, 3. We will now analyze these 
equations for the case of the electroweak phase transition. 

Relevant Scales and Parameters 

To proceed further, we must introduce the relevant physical quantities and 
identify the scales and parameters which are the most useful to the subsequent 
analysis. In the symmetric high temperature phase, we have 



3aT^ 



(4) 



In the symmetry broken phase, we have 



l + L 



T'2 ] 



eh = 3aT^ 



1 - 




(5) 



Here, 



*7r^ 
a = q — 

^ 90 



where g* 107 is the number of particle degrees of freedom. The quantity 
is the critical temperature at which the low temperature and high temperature 

phase have equal free energy, that is, the phase coexistence temperature. We have 
introduced the dimensionless parameter L which measures the fraction of energy 
released in latent heat at the phase transition 



2ar4 



3 eq(re)-eh(re) 

2 eJT,) 



(6) 



L is also a measure of the efficiency of the plasma to absorb the latent heat released 
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at the interface; on this ground, we expect that 



T 



(7) 



The important feature of the electroweak case is the intrinsic weakness of the first 
order phase transition due to the smallness of the parameter L; typically L < 0.01. 
In particular, we can analytically solve for the structure of the deflagration and 
detonation bubbles. 

A second important parameter is the relative velocity of the front with respect 
to the plasma, Vb, which can be written in the form 

IHiT^) = (8) 

where V{T^ = L (1 — Tq/Tc). The function £ is a slowly varying function of 
temperature, which depends on the details such as mean free paths, numbers of 
particles relevant at the temperature of the transition, etc. For our purposes, we 
only need to know the form of V{T), and in addition, we need to know that 
computations of the numerical factors in the above equation, for values of 
relevant for cosmology, give for the EW case ^Vj, ~ 1 with velocities typically in the 
range 0.1 < i'b < 0.9. This range reflects the dependence on the parameters of the 
standard model as well as the uncertainties which result in its determination [23,24]. 
For the QCD case smaller values, Vf, ~ 0.04, are indicated [25]. Unfortunately, more 
accuracy is at present lacking for determining the burning front velocity. 

There are three scales to be considered in the stabihty analysis of the interface 
to follow, and are denoted d^,, dc and i. 

• (ib is a time scale characteristic of the dynamics driving the motion of the 
interface, defined as 

In the electroweak shown in Section 4, df, can be expressed in terms of 
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the thickness of the interface 6 and a parameter e {0 < e < 1) proportional to the 
amount of supercoohng in the system, as (i^ ~ 5/s. 

• dc is the scale over which the surface tension a just balances the difference 
of pressure across the interface, that is, 

dc = . (10) 

Pq -Ph 

Its critical role in the stabihty of the interface was first stressed by Landau [27]. 
In Section 4, it is shown that dc/di^ <^ 1 for the electroweak case, an inequality 
that we will use later. 

• i is the mean free path of the particles which contribute to the heat and 
momentum transfer at the interface. Its effect for the electroweak transition is 
measured by the ratio i/5 which weights the damping of the wall motion in the 
plasma. Because this ratio is typically less than 1, i doesn't appear in the following 
discussions. This is not the case at the QCD scale where this ratio may be larger 
and can play a critical role in the stability of the interface [32]. These diffusive 
effects are not considered here. 

Having introduced the parameters L, i>b and the scales di,, dc relevant to the 
problem and having set their relative magnitudes, we may now proceed. We first 
turn to a detailed study of the kinematics of deflagrations and detonations relevant 
to the electroweak case. In the case of the QCD phase transition, the shapes of the 
bubbles are much the same as is the case for the electroweak phase transition. The 
essential difference is that the latent heat of the transition is larger, L ~ 0.5 — 1, and 
therefore the discontinuities of energy densities are much larger. Since deflagration 
and detonation bubbles have been much discussed in the literature for QCD, we 
shall not repeat this analysis here [20]. 

Kinematics of Deflagrations 



We first turn to the case of defiagrations. Using this parameterization of the 
pressure and energy density, we can now solve for the velocity of the shock front. 
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fgh, the velocity of fluid in between the shock front and the burning front, f flow 5 
and the velocity of the deflagration (burning) front denoted here by v^^f. We 
deflne (see Fig. 3) 

rp2 ^ rp2 ^ rp2 



Tl 



^0 ^ 7-'2 ' ^1 ~ 7-'2 ' ^2 — rp2 



with T2 the temperature to the right hand side of the shock front and the temper- 
ature between the shock front and the burning front is Ti. In the frame where the 
fluid is at rest inside the bubble, the results are 



3ix -j^ 



3(^2 + 2,xl) ' 



^(3x2 + x2)(x2 + 3x2) 



x\ + L(l - xq 



L{l + XQ 



[Zxl + - L(l + xq)] [2>xl + x2 + L(l - xq)] ' 



Xr 



x\ + L(l - xq)] [?,xl + xl + L{1 - xq)] 



[3(x2 - xl) - L(l + xq)] [2,xl + xj- L(l + xq)] 



(12) 



When solving the above equations, we must require that the flux associated 
with the entropy current 



where s is the entropy density, be increased across the surfaces of discontinuity. 
These conditions are the same as 



Xl > X2 

(xQ — xi)^ + L{xq + xi+xqxi — —^xl) — - Xq)>Q 



(13) 



When solving the above equations, we are given the parameters L by knowing 
the latent heat at the phase transition, X2 — 7f by knowing the amount of 
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supercooling from the nucleation calculation, and f^ef by solving the microscopic 
equations for the burning front velocity. 

For the weak deflagrations typical of the electroweak phase transition, we can 
solve the equations explicitly. We have T2 < 3c and assume that Tq — Ti fa 
T2—T\<^Tc. The regions where there are allowed solutions for the above equations 
are then shown in Fig. 4. The linear boundaries of the regions, given in the figure, 
correspond to "^h = "^q = 0) "^h = ''^q = 1) ^"^.L = [22]. The temperature inside 
the bubble is 

T0-T2 L VSVdef 



T2 6 l + x/3^;def' 



(14) 



The temperature in the region between the shock front and the deflagration front 
is 

T1-T2 L ^/ZvA^i 



T2 6 l-3vl 



(15) 



def 



Notice that Ti > Tq > T2 and that (Fig. 4) Tq < Tc- The flow velocity of the 
matter between the deflagration front and the shock front is 

^flow = (16) 
~ "^^def 

There is a singularity in these equations when v^^^ — > 1/3, which is an artifact of 
our approximate solution. When the burning velocity is very close to that of the 
sound velocity in the supercooled phase, then the temperatures are not so close to 
one another. 

Since there is no reason for the burning velocity to be close to the sound 
velocity, the solution above is sufficient for most purposes. For the weak first order 
electroweak phase transition, the temperature is only slightly changed in the fiow 
region and the interior of the bubble wall. The velocity of the fluid is of order v^^^ 
which becomes of order 1 only for deflagration velocities close to that of the sound 
velocity. 
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Kinematics of Detonations 

In the case of detonation bubbles, we have one surface of discontinuity corre- 
sponding to the burning front. Behind the burning front, we must have a similarity 
rarefaction solution of the hydrodynamic equations. We assume that in the region 
of the similarity rarefaction wave, the equation of state of the fluid is close to that 
of an ideal gas. In this case, we can solve the hydrodynamic equations to find 

Q = y-ys 

and 

e = e{ys) e^^^"^^)/^. (17) 

We can now solve for the velocities and energy densities in the various regions 
subject to the constraint of positive entropy generation. In the region of small 
temperature changes we find that 

ei - £2 _ 2 L 

and 

""^""^ 2 37;2 -r ^^^^ 
The energy density far in the interior of the bubble is given by 

^ _ _2 (1 - v^)I 
3 34,-1 

We therefore see that in the energy density in the interior of the bubble is decreased 
relative to that in the outer region. The matter density in the fiow region is 
compressed. Again the changes are small unless the burning front velocity is close 
to that of the sound velocity. In this case our analysis must be modified. 
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3. Stability of Relativistic Planar Combustion Fronts 



In this section we study the stability of planar relativistic combustion fronts. 
The problem has been studied in the non-relativistic case by Landau [27] and 
by Link [26] for non-relativistic velocities but for a relativistic equation of state. 
The result is basically as follows. Small perturbations of deflagration fronts are 
stabilized by the interface tension a while large ones grow exponentially with an 
initial growth time 

^ f b A - Ac 

where Ac is the critical length scale (10) introduced in the previous section 

Ac = 4 . (22) 

We would hke to make two modifications to these earher works. First, in order 
to treat properly the critical role played by the speed of sound Cg ^ l/'s/S, the 
hydrodynamic equations will be written in a completely relativistic form. We can 
thus study the stability for fast deflagrations, w^, Cs and also the stability of 
detonations, > v^^. v^ > Cs. Secondly, we have now available several microscopic 
estimates of the interface velocity v\^, at least for the case — fq <^ [23-25]. In 
particular, this velocity will depend on the temperature Tq ^ < Tc. We shall 
see that if this dependence is sufficiently strong - as seems to be the case for the 
electroweak theory - the conclusions on the stability are completely changed: the 
dependence of v^, on Tq actually stabilizes large scale perturbations of a deflagration 
burning front. 

The stability analysis proceeds as follows. First the leading hydrodynamic so- 
lution, a sharp planar combustion front separating a region characterized by Tq, Wq 
from a region characterized by T^, is described [20]. Secondly, the relativistic hy- 
drodynamic equations are linearized separately in the two regions around a leading 
T — constant, v — constant solution (see also Ref. [28]) and the hnear equations 
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so obtained are solved for the Fourier components of the perturbations of T and v. 
The subtle aspect of the analysis is then the joining of the solutions in the two 
regions across a perturbed interface. The outcome of the stability analysis depends 
on whether or not the connection equations have solutions with — —iuj > such 
solutions grow like exp(ni) and constitute an instability. 

Deflagrations and Detonations 

Here we review the results of Section 2 and formulate them in a form more 
suitable to our purpose. The leading solution is defined as follows. In the rest 
frame of the combustion front (assumed to be in the x, y- plane) , the velocities are 
parallel to the z-axis and the q(h)-phase lies in the region z < (2; > 0) (see 
Fig. 1). We assume that the chemical potential fi = and that the equations of 
state p = p(T) are known. The configuration is thus entirely determined by the 
four quantities Tq, T^, Vq, v\^. These are constrained by the continuity of energy 
and momentum fluxes: 

(23) 

where w = p + e is the enthalpy of the fluid, and of entropy increase 

Sh^hVh > SqTq^^q- (24) 

Thus two of the four quantities remain unspecified by these general relations. 

Information on the further two required relations is obtained from nucleation 

calculations and from microscopic calculations of the wall velocity. The nucleation 
calculations give some temperature Tq < of the q-phase, into which the bubbles 
of the h-phase have to expand. However, as described at length in Section 2, 
boundary condition effects (matter at rest before and after the transition) may 
require the q-matter in front of the interface to be shocked, i.e., have a temperature 
higher than the nucleation temperature. Similarly, a microscopic calculation of 
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gives f h if boundary conditions demand the h-matter behind the front to be at 
rest or if the q-matter in front of the front is at rest. This distinction is not 
important if Vh — 'fq -C Vh- This is the case for the EW transition, as we note from 
the constancy of [Eq. (23)] and the definition of L [Eq. (6)], that 

"q MTc) 2 - ^ ' 

The reason why this velocity difference is so small in the electroweak phase tran- 
sition is because the transition is rather weakly first order, most of the inertia of 
the matter is carried by light mass degrees of freedom which are not much affected 
by the transition and even the W,Z and Higgs bosons gain a mass which is small 
compared to a typical kinetic energy E ~ ST. 

Linearized Hydrodynamic Equations 

We shall now consider either the q- or the h- region and linearize the hy- 
drodynamic equations d^Q'^'^ = 0, Q^'^ = wu^u'^ — pg^'^, = (7,7V), around 
a solution T = constant, v = (0,0, = constant. We shall use as variables 
Sp — Sw/{1 + l/c1), Sv and the transverse 4- velocity variation w^. We assume 
that the sound velocity is the same in the q- and h- phases. Instead of v it is often 
convenient to use the flow rapidity 9 related to it by = tanh^, S9 — ^"^Sv. Thus 
we write 



p^PQ + 5p 

(26) 

u^'^u'^^ Su^" = (cosh 9, 0, 0, sinh 9) + (sinh 9 59, 0, u^, cosh 9 59) 



The components of the zeroth- and first-order variations of the energy momentum 
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tensor are then 



eoo 


= w cosh 6 —p 


5eoo 


= [(1 + cj'^) cosh^ e - l]6p + w sinh 2^ 5^ 


00. 


= ^■u;sinh26' 




= (1 + c^^) 1 sinh 2^ 5p + w cosh 29 69 


0O± 


= 




— w cosh 9 


QZZ 


= w sinh 6 + p 




= [(1 + cj^) sinh^ e + + w sinh 2^ 5^ 


QZl- 


= 




= w sinh ^ u"*" 




= p 


6e^^ 


= 5p 



(27) 



The subscript is omitted from wq and vq whenever no confusion arises. 

The hnearized equations d^,5Q^^ — can then be written in a particularly 
symmetric form by using light cone derivatives d± — e^^(do ± dz). 

i(c;2 + l)d+5p + i(c;2 - l)d-5p + wd+59 + wdx_u^ = 0, 

i(c72 - l)d+5p + i(c72 + l)d_5p - wd-69 + wd±u^ = 0, (28) 

d^6p + lw{d+ + d-)u^ = 0, 

Similarly, the first order variation of the entropy conservation equation d^Ss^ = 
(valid in the bulk fluid, far from the wall) becomes 

^cj^(d+ + d-)5p + ^w(d+ - d-)59 + wdx_u^ = 0. (29) 



Fourier Solution of the Linearized Hydrodynamic Equations 
We shall now search for solutions of Eqs. (28) of the form 

69 = Ae-'''-'^, u^ = Be-'^-^, 6p = Ce''''-'^. (30) 

Since we are interested in solutions behaving as e^* but which do not blow up at 
large distances, it is convenient to introduce Q = —icu and q = ikz and write the 
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exponent in the form 

—ik-x — Q,t + qz + ikx'^. (31) 

The insertion of Eq. (30) into Eq. (28) gives a hnear homogeneous system of equa- 
tions for the Fourier components A, B and C. In general there is a dispersion 
relation 

^(n + vg)2-(g + vn)2 + ^k2 = (32) 

relating the components of the wavevector, and one can express 5p and u-^ in terms 
ol 59. We find two types of solutions: 

dp = -w 69, u = —— 59. (33) 

v\l-\-q viz + 57 

and 

5p = 0, iku^ = ^59. (34) 
7 

The latter is a special solution which describes an incompressible velocity pertur- 
bation (V ■ 5v = 0) moving with the fluid, i.e., it is only a function of z — vt,x-^; 
it is characterized by 

n = -vq. (35) 

The dispersion relation, Eq. (32), deserves further attention; on the {fl,q) plane, 
it appears as the hyperbolae 



q± = z-^ 



with the asymptotes 



Q = q, ^ = -— q. (37 

1 — VCs 1 + CgV 



It is plotted in Figs. 5. From the asymptotes one sees that for v < Cg both 
branches of the curve extend over all values of Jl. However, for v > Cg, one finds 
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q > (actually q > ^/l — v'^c^'-fk) only for 1] < and vice versa. More generally, 
it is a simple exercise to show that the dispersion relation (32) has no solution 
with f > Cg > 0, Re g > and Re ^2 > 0. This fact is of great importance for it 
establishes at once the stability of a detonation front. Indeed, in Section 2, we have 
learned that a detonation front is characterized by > Cs,v\^. On the other hand, 
the perturbation behaves like exp(2;gq); as the q-phase lies in the region z < 0* 
this requires Re > 0. Hence, there is no possible unstable mode (Re > ) 
in a detonation front. In the remainder of our discussions, we will assume that we 
are dealing with deflagration burning fronts only. 

Connecting the Solutions 

We now return to the two-phase situation depicted in the beginning of this 
section and study how its small oscillations behave if the interface is perturbed 
from ^ = to 

z^Cit,x^)^De^^+'^'^^, (38) 

where the amplitude D of the perturbation satisfies 

5<D<l/k, (39) 

The condition that D ^ 5 is just the requirement that we are looking for macro- 
scopic fluctuations in the shape of the wall, that is their characteristic size scale is 
much larger than the thickness of the wall S. The condition that D <^ l/k,l/Q, 
is the requirement that we are looking at small fluctuations. The solutions in the 
q-phase [z < Q and in the h-phase {z > Q can be immediately written down. 
Omitting a common factor of e^*+*'^^ the solution in the q-phase is 

se^ = A, 

K = —F^ A, 

^ + gq 7q 



* Physically imposed after the choice Vq>0. 
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where gq > (to make the perturbation vanish as z ^ —oo). This requirement 
and the fact that we are searching for 1] > formally imply that one cannot include 
the special solution (34) with the dispersion relation Q = —vq in this region, as 
already emphasized by Landau [27]. However, in the h-region we must have < 
and the special solution can (and must) be included: 



5eh^B + C, 

v^il + gh 
1 1 



Uy^ — ik 



B 



v^VL + q^ 7h k2't;h7h 



-C 



(41) 



To see the physics of the special solution, one may compute 



(V X 5vh)* 



k 



7h I ^h^ + ?h 



B + 



1 - 



kvh/ 



c 



in the q-phase no C term appears. As a consequence, in the non-relativistic limit 
V <^ 1 this vanishes in the q-phase and is nonzero in the h-phase: linear per- 
turbations of the interface generate vortices which, through the special solution, 
propagate in the h-phase. 

We thus have four unknowns. A, B, C and D and need four equations to 
determine them (or their ratios). Three of the four conditions are obtained by 
expressing conservation of energy-momentum across the interface d^Q^^ — for 
= ±, 0, z. Here, ©^'^ is the total stress-energy tensor of the system including 
the interface 



H 



-z + c(t, x^)\ + e{;"^ [z - at, x^)] + e^" [<Pk [z - at, x^)] } 



(42) 

The two first terms are straightforward, with Q^'^ and Q'^'^ constant. The third 
term describes the stress-energy arising from distorting the kink solution of the 
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interface, assumed here to be infinitely thin. We have 



with the kink solution satisfying 



Calculating to first order in ( one finds that 

fa{z)-Ca'{z) 

-a{z) + Ca\z) 










-a{z)dxC 



where the notation (for the surface tension) 



a{z)dQC \ 

-(T{z)d^.C, 

-a{z)+Ca\z) -a{z)dyC 

-a{z)dyC / 



a{z)=[M^)]'^a6{z-C), a^j 



dza{z), 



has been introduced. 



(43) 



(44) 



(45) 



If we now compute d^Q^^ from Eq. (42), the result is proportional to 5{z — Q. 
For the two first terms this arises as the derivative of the ^-function. For the kink 
term, Eq. (44) gives d^Qj^ — a{z) S'^^, which again, within the approximation 
(45), is proportional to 5{z — Q. A quantity proportional to a 5 function can only 
vanish if the coefficient of the S function vanishes. We thus obtain, at the wall: 



(46) 



(47) 



The first equations reproduce the leading conservation conditions (23).* The others 
give us three linear relations between the coefficients A, B, C, and D. 



■k In general, there is a term on the RHS of Eq. (46) involving the extrinsic curvature of the 
wall. 
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A fourth equation can be obtained by expressing conservation of entropy across 
the wall, 

Tq^qSq = Th^h^h + J Jsource, (48) 
wall 

where the last term on the right-hand side is a measure of the flux of entropy 
gained by the plasma as it crosses the wall. Prom the work of Refs. [23-25], it has 
been estabhshed that entropy generation reflects directly the existence of a velocity 
dependent force which damps the wall motion to a terminal velocity. This velocity 
was discussed in the second section, 

= ^ . (49) 

This equation of motion has to be modified to include the effect of the curvature 
of the interface as well as of its acceleration, namely, 

a {dl - 9i) c = -F(rq) + . (50) 

Equation (49) yields a unique value for the entropy source in (48) and, consequently, 
can be alternatively used to generate the needed fourth relation. After linearizing 
Eq. (50), we obtain 

v^d^ {di - 9i) C = {Se^ - 7q9oC) - V (51) 

where db — (^/^{Tq), a scale we have already introduced in Eq. (9), and where one 
has defined 

^ V ^dT^J ^q <^^q 

'dTjv^n + q^ ^ ^ 

The latter approximation will be justified later. 
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Equations (47) and (51) form a complete set of homogeneous linear equations 




(53) 



(1 + c, ^hlv^Sp^ + ti;q7q(l + vl){Se^ - doC) 

= (1 + cj^hiv^6pu + w;h7h(l + ^h)(^^h - ^oO- 



(54) 



= [(1 + c-^h^'^l + 1] 5ph + 2^^^h7h^h'^A• 



(55) 



7q^0C + (^o' - ^i)C = (1 - ^) 5^q- 



(56) 



Eq. (53) relates the changes of transverse velocity while Eqs. (54) and (55) describe 
the conservation of the flow of energy and momentum respectively. If rj were zero, 
Eq. (56) tells us that the velocity of the interface relative to the q-phase is only 
perturbed by surface tension effects; this is akin to the boundary condition initially 



on the temperature, this is no longer true and the rj term is needed. Furthermore, 
Eq. (56) provides a direct physical interpretation of the stability of a deflagration 
that we will uncover in the next section. We now proceed to the resolution of this 
system of equations. 

Solution for Small Velocities 

The case of general velocities results in rather lengthy equations and we shall 
first discuss the solution for small velocities. In this case, the equations are much 
simpler and their interpretation is more transparent. This case is defined by 

★ In his analysis, Landau further ignored the effects of the surface tension. This is reasonable 
for most of the applications. 



used by Landau* However, due to the potential dependence of the interface velocity 
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Vn,v^ <^ Cs for which, the integrabihty condition (32) simply imphes 



-k, 



\(jj\ ~ vk. 



(57) 



As a result the solutions in the q- and h-phases simplify to 



5Vr, = A 



5pq = -Wq ( - 



— iA, 



Vq I A 



Svi, = B + C 



-t{B + ^C 



(58) 



After inserting these expressions into Eqs. (53)-(56) and using (23) and (25) as 
well as the definitions of dc, we obtain 



1 JL 



-V, 



n 



LQ. 
q 2 k 



\-{l-V) {^+v^d^k)J 



B 
C 
\DJ 



(59) 



A non trivial solution exists if and only if the determinant of the matrix in Eq. (59) 
vanishes, that is. 



(60) 



or 



1 + 7]- + d],k 



with 



Ac = kz'^ ^dc ( 1 + 



v„ v„ L 

1 + ^ + ^7]- 

^^h ^^h 2^ 



2 



(61) 



(62) 



1-7] I dc 

defined as the critical wavelength. The first solution (60) is spurious and has no 
physical significance. It leads to A — D — B + C — and from Eq. (58) the 
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entire solution vanishes. In general, it lies at f2 = ^^^v\^k, and the solution vanishes 
likewise. The smallness of insures that the term quadratic in Q is small and the 
solution is 



For our purpose, Eq. (63) adequately generahzes the original Landau formula (21). 
Accordingly, the condition of stabihty becomes 



If 77 < 1 the system is unstable for perturbations of scale A larger than Ag. This 
is the case if the wall velocity depends sufficiently weakly on Tq. However, if 
this dependence is strong enough to make rj > 1, large scale perturbations become 
stable. 

We can now comment on two numerical approximations made earlier. Firstly, 
it is obvious from Eq. (63) that fl is of the order of Lv^k so that indeed the 
relevant range of Q is < vk. Secondly, including the dependence of rj in Eq. (61) 
would amount to q — > fq + il/k. However, the additional (positive) term would 
not contribute to the Q-independent terms and thus is irrelevant for the stability 
consideration. 

We now turn to the interpretation of the results (62)-(64). First, we note that 
the critical wavelength Ac [Eq. (62)] receives a correction of order L ^ in respect 
to Landau's value dc [Eq. (22)]; this term originates from the RHS of Eq. (50). In 
a typical macroscopic situation, the latent heat released is a significant fraction of 
the total energy available (L ^ 1) and this correction is neghgible. However, in 
the extremely weak phenomena considered here (L < 1), this correction is large 
and actually controls the value of Ac. 

• If ?7 <^ 1, Ac ^ (4/-^) o?b and is significantly larger than dc because of the 
smallness of L and of the ratio dc/di,. 



n v^{l-r]){l-k/kc) 
k ~ 4 1 + 77 L/2 + 4k 



(63) 



(l-r;)(A-Ac) <0. 



(64) 
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• If T] ^ 1, Ac — —{l/r]) {4:/L) (ib and is driven negative; in which case, 
according to Eq. (64), perturbations at all scales are stable. 

The physics of this stabihty restoration depending on the size of 77 is as follows. 
For the sake of simplicity, let us ignore (with Landau) the RHS of Eq. (51) which 
doesn't play any role in the argument. In the limit 77 — > 0, Eq. (51) reproduces 
Landau's boundary condition. In such a case, the velocity of the perturbation 
9oC follows the variations in velocity of the fluid. This drives perturbations larger 
and larger unless the surface energy d'^C ~ C/^^ is significant enough to prevent 
their growth; this occurs for A sufficiently small. In the present case, the interface 
velocity [Eq. (49)] is proportional to V{T^), an extremely sensitive function of 
the temperature of the q-phase; this sensitivity is measured by rj. The growing 
perturbation triggers a local variation in temperature* which affects, in turn, the 
velocity independent pressure V{Tq) acting on the interface. This increase or 
decrease in pressure acts as a restoring force which opposes the growth of the 
perturbation. When this force is large enough to dominate the others (77 > 1), the 
perturbation collapses. This mechanism operates only if the characteristic growth 
time of the perturbation is larger than the time scale which characterizes the wall 
dynamics, that is, only if ildb < 1. This condition is amply satisfied in the cases 
of interests.^ 

At this point we may also compare with the 4th equation used by Link [26] 
which in the present notation is 



where < a < 1 is a phenomenological parameter used to describe the magnitude 
of in terms of — T^. It is evident that in this case the additional term is 
only a small correction to A and will not modify the conclusions about stability. 
Our calculation, on the other hand, is based on an explicit calculation of Fg = 

★ Positive or negative depending on its direction of growth, 
t As <C k and k is restricted by Eq. (39). 



nA^A+—^a 



( 




(65) 
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^q7q^q(^q) ^^i^d indicates that this dependence may be strong enough to actually 
make large scales stable. 

Solution for Arbitrary Velocities 

As mentioned earlier, the speed of sound Cs plays a critical role in the stability 
of a moving front. In order to extend the previous conclusions to velocities ranging 
up to Cg, one has to generalize the analysis of the previous section. It is convenient 
to introduce the dimensionless quantities 

n = (66) 
Th^hk 



(67) 



The calculation of the determinant proceeds exactly as before and we will only 
quote the results here. Neglecting once again the spurious solution — 1, we may 
write the determinantal equation in the form: 

hv{Q)=Af{Q). (68) 

jV and V are two lengthy polynomials in fi, S'q and S\^. "D is a positive quantity 
unless k/kc is excessively large, while A/" contains the critical information on the 
stability behavior we have uncovered in the small velocity limit. A closer inspection 
of Eq. (68) reveals that the absolute value of f2 ( = [M'/'Dl) is bounded from above 
and is never much larger than one. Consequently, we can attempt to solve (68) 
iteratively in the following manner 

= = etc. (69) 

This iteration converges very rapidly. A comparison of the lowest order of this ap- 
proximation with a numerical analysis shows an agreement better than one percent 
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over the whole range of velocities. We then obtain 



^-7h'^h^ (70) 



with 



P(0) 7h^h (Ah + Aq^^ 7q + 7q^q4k (^Ah + ^Aq) 



+ ^ 2 1 ^q^h^h + Ah7q^q[Aq + (2 + Aq)7qVq7h^;h] 



- 7q^^q4k 



(i + yta;) + ^q^h7q'7h'(2 + a2) (1 + j I 

(71) 



Ac = = + ^4 + ^ ) ^ (72) 



and 



1 - ^ L vAh Aq/ i;q7q 

Solution (70), along with (72), suitably generalizes the corresponding expressions 
(63) and (62) for velocities ranging up to Cg. As a first check, the reader may 
verify that, as and approach zero, S^, S\^, A^ and Aq all approach unity and 
Eqs. (70) and (72) reproduce the familiar small velocity results. Furthermore, the 
condition for stability of a disturbance reads 



Af{n = 0) < (73) 



and it is clear that all the conclusions reached previously hold. In particular, rj > 1 
stabilizes the long wavelength modes. 

As — > Cs, the quantity A^ — > and expression (72) indicates that Ac — > oo in 
this limit; in such a case, our solution (70) is no longer appropriate* However, one 

★ Only for v^^/cs > 0.99. 
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may show, using the full secular equation (68), that, in this limit, the determinant 
possesses no positive roots for any value of k. This result matches nicely with the 
stability of detonations. 

4. Applications to the EW and QCD Phase Transitions 

In this section we apply the results of the previous sections to the evolution of 
an expanding bubble of true vacuum produced at the electroweak phase transition. 
The characteristics of the electroweak phase transition are entirely encoded in the 
finite temperature effective potential for the Higgs field 

V{(p, T) = D{T^ - To^)^^ - ET(P^ + ^0^. (74) 



The parameters D and E may be expressed in terms of the weak coupling 
constant aw = 5'^/47r ~ 1/30 as 

D ^ ^g'^, E ^ (75) 

Here we have used mw ~ ttiz ~ nT't- We can also express the Higgs self-coupling 
as 



At = g'^, (76) 



2 
H 

2 

W 



where rnu is the Higgs mass. This is an effective potential valid near the phase 
transition temperature where multiple scalar degrees of freedom may be important, 
and therefore mn may not be the physical Higgs particle mass [4]. In order to 
generate a baryon asymmetry which is not washed out in the broken phase, the 
effective Higgs mass must be mi/<40 GeV. 

In Eq. (74), Tq is the temperature at which the Higgs mass vanishes and 
is the lowest temperature for which the high temperature phase is metastable. 
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To ~ {m}j/m^r) 250 GeV. The temperature at which the phase transition nucle- 
ates is within the range Tq < T < Tg, where 



r2 - r2 _ £;2 



T2 \tD 

^awml, (77) 

~ 2~ 

TT rrijj 
<0.04. 

In terms of these parameters, 

Ph = aT^ - V{T^) , (78) 

where V{T\^) is the value of ^(0, Tp,) at its absolute minimum. We can adequately 
express it as 

V{T^) = eV, (79) 



with 



where e measures the amount of supercooling at the phase transition; it has the 
value 0.25 for m/y = 40 GeV. 

The wall thickness 6 can conveniently be expressed in terms of Vc and the 
surface tension a — J dcf) ^/2V{(f), Tc) as 



wall ^ ^ ^ (81) 



~ — ~ 10- 100 T-^ . 



We can now develop an intuition for the scale db {— (^/^{Tq)) which appears 
naturally in our calculations; from its definition and Eqs. (79) and (81), we can 
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express it as 



4 ^ - ■ (82) 

e 



In the situation of maximal supercooling {e = 1) = 5 and it increases as the 
phase transition becomes more weakly first order [e 0). 

According to Eqs. (5) and (78), the relative latent heat is 

- _ 8 E'^D ( ^ ^2 ^^ 



2a A| 

5 „2 (83) 

4 
H 



<0.01 



In the case of 77 > 1, we have asserted that there are no unstable perturbations 
on the basis that perturbations of size larger than the critical wavelength Ac are 
stable and that Ac has the negative value 

Ac^---4- (84) 
^ L 

Prom Eq. (62), this is a trivial assertion assuming L < 1 and the ratio dc/d\, is 
smaller than one. We now prove the latter assumption. Remembering that dc is 
defined as cf/{p^ — Ph)) using Eqs. (5) and (14) as well as (79) and (83), this 
ratio satisfies 

dh Pq - Ph Pq - Ph 

To compute the parameter 77, we only need to know the temperature dependence 
of the interface velocity v-^, which is [23,24] 

Vh^ e = _ . (86) 
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The parameter tj [Eq. (52)] then becomes 



V 




(87) 



e aw rn^' 



As £ behaves parametrically as m^j , we can obtain a lower bound for 77 by 
considering the lowest value for ttih, which, to be conservative, we take to be 
40 GeV, for which e ~ 0.25 [23,24]. In such a case. 



and unstable perturbations can only develop for 



a range clearly not favored by the latest estimates: 0.1 < Vb < 0.9. Prom the 
numerical results above, we find that perturbations on all scales are unable to 
destabilize the shape of a bubble: if produced spherical, it will stay spherical until 
it collides with other bubbles and completes the transition. This conclusion, which 
contradicts a previous analysis [31], is illustrated in Figs. 6a and 6b. 

For the QCD quark-hadron transition, the problem is entirely nonperturbative 
and results are much less clear. L is expected to be larger, in the range 0.5-1. 
In the cosmological context Vb is rather small [25], of the order of 0.04. A fair 
calculation of rj is at present an impossible task. From Eq. (87) we see that unstable 
perturbations can only develop for 



The standard nucleation analysis for the cosmological phase transition [25], using 
values of L and a measured with lattice Monte Carlo techniques [29-30], suggests 




(89) 




(90) 
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1 — T^/Tc ~ 0.001 or an upper limit on of about 0.03. This is very close 
to the estimated value 0.04 of fb and thus no definite statement can be made. In 
terms of the present analysis, QCD, quite interestingly, seems to be a borderline 
case. 

5. Conclusions 

The prospect that baryogenesis took place at the electroweak scale has recently 
been the motivation for the investigations of various aspects of the electroweak 
phase transition (EWPT). In the present work, we have analyzed the structure 
of a propagating front separating the unbroken phase from the broken phase. In 
particular, we have generalized the analysis of Landau [27] of the hydro dynamical 
stability of a moving front under small perturbations. We have found that, due 
to the particular nature of the EWPT, this generalization leads to a conclusion 
opposite to previous works on the subject [31], namely, that the growth of an 
electroweak bubble is stable under small perturbations. 

This surprising result comes about because of the intrinsic weakness of the 
EWPT, which is characterized by the smallness of the ratio of scales dc/cib (~ Oiw) 
and by the smallness of the latent heat relative to the energy flowing across the 
interface L ~ «^(90/7r^(7*), which, in particular, result in a large sensitivity 
rj ~ 1/aw of the velocity of the interface on the temperature of the flowing 
plasma. In general, there is a critical velocity below which instabilities are allowed 
to develop. In the EW case, we found that this critical velocity is parametrized 
as ~ ^/awj^ {rnw / tTi^h) and is bounded from above by 0.07. This value is 
significantly smaller than the values obtained from microscopic calculations [23,24], 
0.1 < Vb < 0.9. 

It is interesting to ask under what conditions a smaller value of could be 
obtained. In the minimal standard model, it would happen [23,24] if the mean free 
paths £ of the gauge bosons and of the top quarks are larger than the thickness 
of the bubble wall, d\^. The velocity can in such a case reach a value of 0.04 
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{nhop ~ 140 GeV) and is slightly smaller as the top quark mass increases. In a 
non-minimal standard model, smaller values might be allowed because of the slight 
dependence of on the microscopic parameters. 

In the QCD phase transition, too little is known to make definite statements 
on the stability of the quark-hadron interface. One must remember that even 

the order of the phase transition is not known at present. It is quite interesting 
that reasonable estimates of the QCD parameters invoked make the cosmological 
quark-hadron phase transition a borderline case of our stability analysis. 

The authors wish to thank the Institute for Theoretical Physics at the Univer- 
sity of California, Santa Barbara, where this work was begun. 
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Figure Captions 

Figure 1: The burning front in the rest frame of the front. The wall is moving 
to the left, the inflowing velocity of the high temperature phase is and the 
outflowing velocity of the low temperature phase is v\^. Deflagrations are when 
< fh and detonations when > v^. 

Figure 2: Structure of a detonation bubble in 1+1-dimensions: (a) The time 
dependence of the x > part of the bubble with the three fronts demanded by 
the boundary conditions. The dotted line shows the path of a fluid element, (b) 
The energy density and flow rapidity of matter as functions of space-time rapidity 
y — tanh.~^{x/t) for a detonation bubble. For Jouguet detonations yi — ys and no 
region of constant velocity appears. The magnitudes of the quantities shown are 
related by Eqs. (17)-(20). 

Figure 3: Structure of a deflagration bubble in l+l-dimensions: (a) The time 
dependence of the x > part of the bubble with the two fronts demanded by the 
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boundary conditions. The dotted line shows the path of a fluid element, (b) The 
energy density and flow rapidity of matter as functions oi x/t for a deflagration 
bubble. The magnitudes of the quantities shown are related by Eqs. (14)~(16). 

Figure 4: Regions allowed by the continuity of energy and momentum fluxes 
and entropy increase for combustion from a high T phase at T = Ti to low T phase 
at r = Tq. The three approximately hnear boundaries follow for Ti — Tc, Tq — Tc <S 
Tc from Eqs. (12) by demanding the velocities to be real and from the second of 
Eqs. (13). The curves fdof = const, are explicitly given by Eqs. (14) and (15). Note 
that detonations are possible only after supercooling by at least L/4. Horizontal 
hatching corresponds to the forbidden region characterized by As < 0, and vertical 
hatching to imaginary velocity. 

Figure 5: A plot of the dispersion relation (36) for (a) v < c^, (b) v > Cg. 

Figure 6: Long wavelength limit of Q/k computed for the electroweak theory 
(m/i = 60 GeV). (a) Landau's results (21) (dashed curve) and our results (63) 
(solid curve) are shown for small velocities v <^ Cg. (b) In the limit v — > Cg, the 
inadequacy of (21) and (63) is illustrated by comparison with our general formula 
(70) (dotted curve). 
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